Fano Effect through Parallel-coupled Double Coulomb Islands 



in 
o 
o 

(N 

> 
O 



on 



i 

c 

o 
o 



> 
o 



in 
o 
-i— > 



i 

C 

O 

o 



Haizhou Lu/'Q Rong Lii, 1 and Bang-fen Zhu 1,2, Q 

1 Center for Advanced Study, Tsinghua University, Beijing 100084, China 
2 Department of Physics, Tsinghua University, Beijing 100084, China 

By means of the non-equilibrium Green function and equation of motion method, the electronic 
transport is theoretically studied through a parallel-coupled double quantum dots(DQD) in the 
presence of the on-dot Coulomb correlation, with an emphasis put on the quantum interference. It 
has been found that in the Coulomb blockage regime, the quantum interference between the bonding 
and antiboding DQD states or that between their Coulomb blockade counterparts may result in the 
Fano resonance in the conductance spectra, and the Fano peak doublet may be observed under 
certain non-equilibrium condition. The possibility of manipulating the Fano lineshape is predicted 
by tuning the dot-lead coupling and magnetic flux threading the ring connecting the dots and 
leads. Similar to the case without Coulomb interaction, the direction of the asymmetric tail of Fano 
lineshape can be flipped by the external field. Most importantly, by tuning the magnetic flux, the 
function of four relevant states can be interchanged, giving rise to the swap effect, which might play 
a key role as a qubit in the quantum computation. 

PACS numbers: 73.23.Hk, 73.63.Kv, 73.40.Gk 



I. INTRODUCTION 

Fano resonance stems from quantum interference be- 
tween resonant and nonresonant processes,- and mani- 
fests itself in spectra as asymmetric lineshape in a large 
variety of experiments. It is known that Fano effect is a 
good probe for the phase coherence for carriers in solids, 
in particular in a quantum dot (QD) systemiSiMiiiiSi 7 . Un- 
like the conventional Fano resonance )S*2ii2iii the Fano 
effect in the QD system has its advantage in that its 
key parameters can be readily tuned. Suppose a discrete 
level inside the QD is broadened by a factor of T due 
to coupling to the continua in leads, the key to realize 
the Fano effect in the conductance spectra is that within 
r centered at the resonance level, the phase of the elec- 
tron in the non-resonant channel changes little^Siii The 
first observation of the Fano lineshape in the QD system 
was reported by Gores et al?'^ in the single-electron- 
transistor experiments. Recently, K. Kobayashi et al. 
studied the magnetically and electrostatically tuned Fano 
effect in a QD embedded in an Aharonov-Bohm(AB) 
ringj££ and A.C. Johnson et al. investigated a tun- 
able Fano interferometer consisting of a QD coupled to a 
one-dimensional channel via tunneling and observed the 
Coulomb-modified Fano resonances. 7 

Recent experimental advances in the parallel-coupled 
double quantum dots (DQD) 14 i 15 i 16 i 17 i 18 , in which two 
coupled QD's via barrier tunneling are embedded into 
opposite arms of a AB ring respectively and also coupled 
to two leads roughly equally (cf FigQJ, have inspired 
a number of theoretical attempts to study the coher- 
ent transport in this systemi^&SliS^l As a controllable 
two-level system, it is appealing for the parallel-coupled 
DQD system to become one of promising candidates for 
the quantum bit in quantum computation based on solid 
state devices^* The entangled quantum states required 
for performing the quantum computation demand a high 
degree of phase coherence in the system^ Being a probe 



to the phase coherence^ the Fano effect in the parallel- 
coupled DQD system is certainly of practical importance, 
if tunable, and especially if the swap effect can be ma- 
nipulated. 

In a real small system like quantum dots, the electron- 
electron (e-e) interaction will influence the transport pro- 
cess. A question naturally arises, whether the e-e inter- 
action in the QD breaks the Fano resonance in the DQD 
system or not? Till now, most theoretical works address- 
ing the Fano effect in DQD ignored this aspect by adopt- 
ing a Fano- Anderson model in which the e-e repulsion is 
entirely absent ^S&SLS^i^ The electron correlations in 
the DQD have been taken into account in coherent elec- 
tronic transport through a DQD connected in a series 
with electrodes jSl however, no Fano resonance has been 
found in this configuration. In this article, we intend to 
investigate the Fano effect in the tunneling-coupled paral- 
lel DQD system in the presence of finite on-dot Coulomb 
interactions, that is, the Fano effect in Parallel-coupled 
Double Coulomb Islands. 

The paper is organized as follows. After introducing 
the two-impurity Anderson model with an inter-dot cou- 
pling term to describe the parallel-coupled DQD system 
in Sec[n] the current through the mesoscopic system is 
formulated in Sec with the Green functions in the cen- 
tral region and Fermi distribution in leads. 28 In small U 
cases, all spin flip processes may be ignored, within the 
Hatree-Fock approximation(HFA), the current formula is 
reduced to a Landauer-Biittiker one, and the electrons 
with spin a behave like moving in a mean field of elec- 
trons with the opposite spin W. Then numerical results 
for the zero temperature case are presented in Sec II VI 
Owing to the on-dot repulsion and inter-dot coupling, 
there are four quantum states in the DQD system rele- 
vant, centered approximately at Eq — t c , £o+t c , £q — t c + U 
and £o + t c + U, corresponding to the bonding and anti- 
bonding DQD states and their Coulomb blockade coun- 
terparts, respectively. It has been found that the Fano 
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resonance appears as a result of the quantum interference 
not only between the DQD states at e — *c an d e + t C7 
but also between their Coulomb blockade counterparts. 
It has also been found that the direction of the asymmet- 
ric tail of Fano lineshape can be flipped by the external 
field. Most importantly, by tuning the total magnetic 
flux through the AB ring, the swap effect between four 
resonance peaks in the conductance spectra is predicted, 
which might be of potential application as a type of C- 
Not gate in the quantum computation. Finally, a brief 
summary is drawn and presented. 



II. PHYSICAL MODEL 



The total Hamiltonian is expressed as 



H = H, 
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FIG. 1: Schematic setup of a tunneling-coupled parallel DQD 
system coupled to two reservoirs. 

We start with the two-impurity Anderson model for 
the parallel-coupled DQD (c/Fig0 as 

Hdqd = >J Siydlpdig + Uid^d^dlydii 

~ tcY,( d Ld2* + h.c), (2) 



where for the sake of simplicity the tunneling matrix 
element V a i is assumed to be independent of k, and 
the phase shift due to the total magnetic flux thread- 
ing into the AB ring, $, is assumed to distribute evenly 
among 4 sections of the DQD-AB ring, namely Vli — 



\V L1 \e l ~,V 1 



L'l 



\V L 2\e^,Vj 



|Vm|e l 4, and V R2 



\Vn2\e lT ■ Here <fi — 27r$/<f>0! i n which the flux quan- 
tum $0 = he/ e. In the following calculation, we define 
the line-width matrix as Tf- = J^k VaiV*j2n5(e — Ek a ) 
(a = L, R) and T = T L + T R . According to FigQ the 
line-width matrices in the QD representation read 
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where Tf is short for 
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III. CURRENT FORMULA 

According to Meir and Wingreen^ the general for- 
mula for current through a mesoscopic region between 
noninteracting leads is given by 



dwTr{(r L - T /? )G < 



+ [f L (uj)T L - f R (oj)T R ](G r (u J ) - G»)},(6) 

where, /l(_r)(w) is the Fermi distribution function on the 
left(right) leads, T L< - R ^ has been given by Eq.©, G r , 
G° and G < are the retarded, advanced and lesser Green 
functions in the DQD region, respectively. 



where d\ a idia) represents the creation (annihilation) op- 
erator for the discrete state with the energy Ei a and spin 
a ( a =|, I ) in dot i ( i = 1, 2), which are coupled each 
other via tunneling t Cl and the on-dot Coulomb repulsion 
is described by the second term in Eq.Q. 

The Hieads in Eq.Q represents the noninteracting elec- 
tron gas in the left(L) and right (R) leads, 



(3) 



where c kaa (ckaa) is the creation (annihilation) operator 
for a continuous state in the lead a (a = L, R) with 
energy Ska and spin a. 

The Ht in Eq.QJ represents the tunneling coupling 
between the QD and lead electrons, 



H 
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A. Retarded and Advanced Green Functions 

The retarded Green function is defined as 

G^ ja {t) = {{dUt)\d] a )Y = -i0(t)({d itT (t),dl}). (7) 

Writing the equation of motion for the retarded Green 
function in Fourier space {{di^d^)) 1 ^^ one arrives at 

(w - e ia + l -r u )((d i<T \d a ))i + (t c + l -r u )((d la \dU)i 



Sij + Uiddi^niald^}^. 



(8) 



where, i = 2 if i = 1, and vice versa. During the deriva- 
tion, we have calculated 

// ut \\r Vaj((djcr\dl a )) , ~j 
{{Ck<xa\d} a }} r u = 2^ ( 9 ) 



i=i,2 



W - Ska + i0+ 
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and defined the retarded self energy originated from the 
dot-lead couplings as 



A: , ft 



(10) 



where + represents an infinitesimal. In the wide-band 
limit. 



and 



/OO J 



(17) 



The derivation of advanced Green function G a follows 
the same procedure, and its expression is the Hermite 
conjugate of the retarded Green function. 



yr 



-_(T L + r R ) = ~-r 



(11) 



The equation of motion for the last term in Eq.JSJ), 
((dianiald]^))^, reads 

(w - e ia - U i ){{d i(T n i ff\d] (T ))l } 
= Mn*r> - t c {{d lrT n m \d) a ))l - t c {{d ta dU^\d] a ))l 

+ UidiadLdi^d]^ +J2 V ci(( c kaan i ^\dt a )) r 0J 

k.a 



k,Oi 

Y, v aM a *dMdi)Z- 



(12) 



k.a 



To truncate the set of equations for the retarded Green 
function, we adopt the HFA to the higher-order Green 
functions on the right side of Eq. (|12(l and have 

m^idU): s (n^adjjdUz, 



((d ia did~\di))i - (4^>((^i<»^ 

{(ck a an lW \d] a ))l =i (n t w)((c ka *\d] a ))l, 
({d) m Ck«Md\))l - (dlc kaW )((d ia \dU)l. (13) 



Thus obtained retarded Green function is expressed in a 
compact form as 



G» = [l-g r (w)El- 1 g r («) 



(14) 



in which S r is given by Eqs. iJSJ and l|ll|) . and g r (w) is 
the Green function for the isolated DQD. It is convenient 
to express the inverse of g r (w) as 



and 



(u - e icr )(io - e ia - - Uj) 

to - E ia -Ui + Ui(n l7 j) 

| Ujt c [{dU~) - {cLdit)] 



(15) 



Here, the expectation values of (n^) and (d\—dj-) can be 
calculated self-consistently by taking advantage of the 
definition of the lesser Green function 



(riiw) = -i \ —G ilTi& {uj), 



(16) 



B. Lesser Green Function 

The lesser Green function is defined by 

Gt, ja {t) = ((d ia (t)\dl))< = i(d] a d ta {t))- (18) 

The equation of motion for the lesser Green function in 
the Fourier space follows as 



{u-e ia ){{di.\d\ a ))<+t c {{d.y ja ))< 



V al {{c kaa \d]a))t+U l {{d l an l ^\d) a ))<. (19) 



t \\< 



k . a 



The first term on the right hand side of Ea. l|19() can be 
achieved with the analytic continuation^ as follows 

((c ka «\dl))< = Ki 



1=1,2 
vt \\< 



x[ 5 LM((^M] CT })<+^ a H((d iCT |dT a ))2] ) (20) 

where £fl„ is the Green function for the noninteracting 
electrons in the leads, 



#L( W ) 



1 



W — Eke 



£0+' 



dkai^) = if a (u)27r5(u - e ka ). 



(21) 



The equation of motion for the last term on the right 
hand side of Eq. 1|19|) has almost the same structure with 
Eci. (|12|) except for term dij(riia). Under the Hartree- 
Fock approximation, by using the expressions similar to 
Eci. (|13|) . we eventually arrive at 



(uj - E icr - Ui)((dia-n iw \d'.)) 



t \\< 



tc(n iw )((d la \dU}< - t c (dUia)((dia\d J j<7 )) 



t \\< 



+ tMM((d t M a ))< 



1=1,2 



+i(fLTt l +f R T*)((d la \d] a ))Z]. (22) 



Inserting Eqs. (|20|l and into Eq. (jT§|l . the expression 
of G < can be simply cast into 



G< = G r S < G t , 



(23) 



where £< = i{fiX L + JrT r ), indicating that, with 
the second-order HFA, the self energies due to the dot- 
lead coupling are separated from those resulting from 
Coulomb interaction. 



4 



C. Current Formula 

Generally, G' G Q = G r (S' r - S a )G a . With the 
discussion above, the Eq. JJJJ of current is simplified to2& 

J = J2lJ du ) [f L (u)-f R (w)]Tr[G a (w)r R G r (u)T L ], 

(24) 

where f L{R) = l/{exp[(w - /i T ^-)/k B T] + 1}, and V is 
the applied bias voltage. This current expression reduces 
to a usual Landauer-Biittiker formula for the noninter- 
acting case, implying that, with up to the second order 
of the HFA, the system can be effectively described by a 
single-particle picture. In this context, the effect of elec- 
trons with spin a on the motion of the electron with spin 
a behaves like a background, and the coherent tunneling 
process takes place between electrons with the same spin. 



I T)i| Ik I Dil T>2, | U)i|0) 2 , |0)i| T|> 2 . The ground 
and excited two-electron-states are determined by di- 
rectly diagonalizing the matrix representation of Hdqd, 
spanned in the Hilbert space by these 6 states. Thus the 
two-electron ground state is associated with an energy at 
2e + ±(U- ^/U 2 + 16*2) w 2e -At 2 jU ~ 2e - Since the 
intra-dot Coulomb interaction produces effective charg- 
ing energy on the bonding and antibonding states, it is 
expected for the two-electron ground state that the elec- 
trons tend to distribute themselves evenly throughout the 
DQD structure to avoid the charging energy. 

Because an extra charging energy has to be consumed 
to add the third and fourth electron into the system, then 
the ground states energy is approximately 3eo + U — t c 
and 4eo + 2U, respectively. 



IV. NUMERICAL RESULTS AT ZERO 
TEMPERATURE 

The parameters for our model calculation are taken as 
follows. For simplicity, we assume £\ a — e^o = £a = 0. 
Throughout this work, we consider the case where the on- 
dot charging energy U is much larger than the interdot 
coupling t c which is of the order of unity, thus we take 
U-y = U% = U = 4. As pointed out in Ref 12(1 in the 
parallel-coupled geometry, if rf = T% and Tf = rf , 
the antibonding state could be decoupled entirely from 
the leads. Our calculation verifies that this situation is 
retained when the on-dot Coulomb repulsion is taken into 
account. To avoid it, in the following we only choose 
two configurations: (1) rf = rf > T% = rf , and (2) 

T^L -nR ^ pi j^R 

1 1 — 1 1 > 1 2 — 1 2 ■ 

It should be pointed out that the numerical results 
presented in this Section are valid only at temperature 
above the Kondo temperature, though are calculated at 
zero temperature, because all the spin flip processes have 
been neglected as mentioned above. 



A. Equilibrium process 
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We are particularly interested in how the states in 
the DQD region are modified by the on-dot Coulomb 
repulsion, and how this modification influences trans- 
port properties. First, let us estimate the eigenstate and 
eigenenergy of the isolate DQD system. 

In the case of only one electron in the DQD system, 
due to the inter-dot coupling, the DQD states are the 
linear combination of the states in two dots, thus the 
formed bonding and antibonding states are associated 
with energy at £q — t c and eq + t c , respectively. Hence, 
the one-electron ground state of the DQD is the bonding 
state with eigenenergy Eq — t c . 

When the DQD contains two electrons, 6 possible 
states in the system include: | O2, | 4)2, 



FIG. 2: (Color online) (a)The local density of states and (b) 
the total particle number in the DQD structure for different 
magnetic phase parameter <j>. A magnetic flux threading into 
the AB-ring swaps the effective couplings of the DQD states 
as well as their Coulomb counterparts to the leads. The pa- 
rameters taken for calculation are: e\ a = Eio = 0, f = 1, 
Ui = U 2 = 4, rf = If = 1 and if = rf = 0.15. Same 
results are obtained when the dot-lead coupling is changed to 

rf = rf = 1 and rf = rf = 0.15. 

When an isolate DQD is connected to two leads, the 
couplings between the DQD levels and two leads re- 
sult in broadening of the discrete levels and forming 
bands. The local density of states (DOS), defined as 
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the imaginary part of the retarded Green function p a = 
— — 2t=i 2 I m( ^T<T i<r I have been calculated for three dif- 
ferent magnetic flux 4>. Fig. Ufa) reveals that the line- 
shape of the local DOS critically depends on magnetic 
flux 4> through the modulation of the effective coupling 
between the DQD states and leads, which can be tuned 
by both the dot-lead coupling strength and the total mag- 
netic flux. Fig.[2a) also shows the band width variation 
versus the magnetic flux: The broadening of the bonding 
DQD state is always accompanied by the shrinking of the 
antibonding DQD state; same for their Coulomb coun- 
terparts, though the lineshape of the DQD states and 
their Coulomb counterparts are somewhat different. It is 
also noted from Fig. |2ta) that, however, unlike without 
Coulomb repulsion case2£ where the total width of the 
bonding and antibonding bands is an invariant because 
the self-energy is solely determined by the DQD-lcads 
coupling; with the on-dot Coulomb correlation into con- 
sideration, an additional self-energy due to Coulomb re- 
pulsion plays a role, then the total band width in general 
depends on the magnetic flux <f> to some extent. 

According to Ea. (|24|l . the differential conductance at 
equilibrium is defined as Q(n) 
temperature reads 
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FIG. 3: The differential conductance for configuration 1 as a 
function of Fermi level. The Fano resonances are marked with 
arrows. The parameters taken are: s\ a = £20- = 0, t c = 0.8, 
rf" = If = 1, rf = If = 0.15, and Ui = U 2 = 4. 



(25) 



where /i is the Fermi level. When Fermi level varies, 
the occupation number of electrons in the DQD re- 
gion is changed correspondingly. As shown in Fig^b), 
the integer number of electrons confined to the DQD 
region occurs approximately at the following energies: 
£0 - t c , e + t c , e + U- t c , e + U + t c M£L 

Figs. and 0] show the differential conductance spec- 
tra as functions of the Fermi level (or equivalently, of the 
dot level Eq) in two different configurations. The peaks 
marked with arrows represent the peaks of the Fano- type, 
whose lineshape is asymmetric compared with the sym- 
metric Lorentzian at the same spectra. Four Lorentzian 
peaks appear in the conductance spectra when <p = n; 
while there are only two Lorentzian and two Fano peaks 
for <p — 0. The spectra for (j) — lix has a mirror symmetry 
with that for </> = 0. 

The peak positions at the DOS and conductance spec- 
tra we obtained are in good agreement with those in 
Refs. l3ll and l2?l in which two dots coupled in series were 
considered; however, no Fano effect was observed in the 
in-series-coupled DQD system. It is then necessary to 
explain why the Fano interference occurs in the parallel- 
coupled DQD system, while absents from the in-series 
DQD system. 

As pointed out in our recent article^ in general, the 
lead states couple to the bonding and antibonding levels 
with different strengths, leading to quite different broad- 
ening. If the band width of the strongly-coupled level 
covers the weakly-coupled band and the phase shift for 
states in the strongly-coupled channel is negligibly small 
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FIG. 4: The differential conductance for configuration 2 as 
a function of Fermi level. The Fano resonances are marked 
with arrows. The parameters taken are the same to the figure 



above except for t c = l,Ti 



1 2 



1, and rf 
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0.1. 



when across the weakly-coupled level, then the Fano in- 
terference will occur, in which the weakly-coupled chan- 
nel acts as a Breit-Wigncr scatter in the resonant tun- 
neling process while the strongly-coupled channel can be 
regarded as a non-resonant one. For example, the trans- 
mission amplitude of an electron resonantly traversing 



6 



through the a weakly-coupled energy level at eo can be 
described 



tR- 



u - eo 



■iT' 



(26) 



which implies a phase shift by tt for the transmission am- 
plitude over an energy range of T around this level, and 
tji is a coefficient. On the other hand, the transmission 
amplitude for states in the strongly-coupled channel can 
be approximated as tMe l ^ N , in which the phase 4>n varies 
little when across the narrow band centered at eo, the in- 
terference between these two channels yields the Fano 
lincshapc around the weakly-coupled level eo as 



So + iT 



+ t N e 



i<t>h 



1 



(27) 



where the detuning e = (u> — eo)/T, and the asymmetric 
factor q = i + tfie~ 1 ^ N /t^ . It should be noted that in the 
present system, the resonant or non-resonant channel is 
not fixed, on the contrary, it could be any one of the four 
channels, as long as the required magnetic flux as well as 
Fermi level are satisfied. 

Notice also that when the dot-lead coupling strength 
is adjusted such that the configuration 1 (FigJSJ is trans- 
formed into the configuration 2 (Fig^J), the tail direction 
of the Fano peak is flipped. The reason behind this be- 
havior is that, compared to configuration 2, an extra flux 
of 7r threads into the loop in the configuration 1^ then 
the Fano lineshape in configuration 1 is just opposite to 
that in configuration 2. Namely, if two channels inter- 
fere with each other constructively in configuration I, 
they will interfere destructively in configuration 2, and 
vice versa. Because of different inter-dot coupling t c , the 
splitting between the bonding and antibonding bands in 
Figs. Eland 0| is different, resulting in different Fano line- 
shapes in these two figures. 

As for the in-series DQD geometry, the bonding and 
antibonding states are coupled to the leads with equal 
strength since one dot has to couple to the left or right 
leads via the other dot. Though the quantum interference 
can occur in this geometry as well, one DQD state can 
not cover or be covered by the spectrum of the other state 
to give rise to the Fano resonance. 



B. Out of Equilibrium 

In the presence of a finite bias voltage V, the differen- 
tial conductance Eq. (|24|l becomes 



a a=L,R 



Tr[G a ^ a )T R G r ^ a )T L 



(28) 



where hl(r) = and/i = {^ L +pb R )/2. Varying /x is 

equivalent to adjusting the DQD energy levels reversely. 
Eq. I|28|) indicates that the left and right leads contribute 
to the differential conductance separately. Thus, the 
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FIG. 5: (Color online) The conductance as functions of the 
average of the left and right Fermi levels, /i, for different bias 
voltages in configuration 1. (a) The resonances at eo —t c and 
£o + t c for <f> — 0. (b) The resonances at £o + U — t c and 
eo + U + t c for <j> = 2tt. 



Fano peak doublet in the conductance spectra can be ex- 
pected, when either of two Fermi levels is aligned with the 
weakly-coupled DQD level under certain circumstance. 
More specifically, when the system is driven by a bias 
V, if eV is comparable with or greater than the width 
of the weakly-coupled level, but less than the width of 
the strongly-coupled level, the original Fano peak in the 
conductance spectra may be split. In Fig0 two exam- 
ples are presented: (a) The conductance spectra around 
Sq — t c and Eq + t c in the absence of magnetic flux, and 
(b) the spectra around Eo + U — t c and Sq + U + t c in 
the presence of a 2ir flux. In both cases, three values 
of eV locate within a region, which is comparable with 
the width of the weakly-coupled states and less than the 
width of the strongly-coupled states. It is obvious from 
the Figure that the splitting of the Fano doublet is pro- 
portional to the applied bias, which, together with the 
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Lorentzian peak, makes a step-like lineshape in the con- 
ductance spectra. 

V. CONCLUSIONS 

In summary, within the Keldysh non-equilibrium 
Green function formalism and by using the equation of 
motion method in which the Hatree-Fock approximation 
is applied to the higher-order Green function, the trans- 
port through the parallel-coupled DQD system has been 
studied with an emphasis put on the effects of the intra- 
dot Coulomb correlation on the Fano interference. We 
have predicted that the Fano effect, as a consequence 
of quantum interference, could survive in the presence 
of a finite Coulomb repulsion. Without loss of general- 
ity, we select four states corresponding to the bonding 
and antibonding DQD states and their Coulomb block- 
ade counterparts as the basis for investigation. The cou- 
pling between each of these four states and the electrodes 
broaden the level into a band with different width. The 



mechanism of the Fano lineshape in conductance spectra 
is explored. It has been also found that the direction of 
the asymmetric tail of Fano lineshape can be flipped by 
tuning the dot-lead couplings. More interestingly, by ap- 
plying a magnetic flux, linear responses of the four states 
could be interchanged, leading to a magnetic flux tun- 
able Fano effect. When a suitable bias voltage is applied, 
the Fano peak doublet and a step-like lineshape may be 
observed as a result of the splitting of equilibrium Fano 
resonance. 
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